Low-energy properties of anisotropic two-dimensional spin-1/2 Heisenberg models in 

staggered magnetic fields 
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We present a systematic study of the anisotropic spin-1/2 Heisenberg model in staggered magnetic 
fields in two dimensions (2D). To mimic real materials, we consider a system of coupled, antiferro- 
magnetic chains, whose interchain interaction can be either ferro- or antiferromagnetic. When the 
staggered field is commensurate with the magnetic interactions, an energy gap opens immediately 
and follows a power law as a function of the applied field, similar to the situation in ID. When 
the field competes with the interactions, a quantum phase transition (QPT) occurs from a gapless, 
magnetically ordered phase at low fields to a gapped, disordered regime. We use a continuous-time 
Monte Carlo method to compute the staggered moment of the ordered phases and the spin gap of 
the disordered phases. We deduce the phase diagrams as functions of the anisotropy ratio and the 
applied field, and calculate the scaling behavior of the models in both quantities. We show that in 
the competitive case, the staggered field acts to maintain a regime of quasi-lD behavior around the 
QPT, and we discuss as a consequence the nature of the crossover from ID to 2D physics. 

PACS numbers: 75.10.Jm, 75.30.Kz, 75.40.Cx, 75.40.Mg 



I. INTRODUCTION 

The spin-1/2 antiferromagnetic Heisenberg model in a 
staggered magnetic field has attracted increasing interest 
as the number of real materials it describes continues to 
rise. While the model in zero field is gapless in all dimen- 
sions, the staggered field lowers symmetries, changes uni- 
versality classes, and induces gapped states and uncon- 
ventional, gapped elementary excitations. Intensive in- 
vestigations into staggered-field effects in quantum mag- 
nets were motivated originally by inelastic neutron scat- 
tering experiments on the quasi- ID spin chain copper 
benzoate, 1 Application of a uniform magnetic field to this 
material caused an unexpected gap to open in the exci- 
tation spectrum. This energy gap follows a power law 
as a function of the applied magnetic field, but its mag- 
nitude depends strongly on the direction of the applied 
magnetic field. 

This unusual experimental finding was soon 
explained 2,3 by the fact that copper benzoate is 
not a perfect chain, but has an alternating crystal 
structure giving rise to a staggered gyromagnetic tensor, 
possibly accompanied by a Dzyaloshiskii-Moriya (DM) 
interaction. The full Hamiltonian in a uniform external 
magnetic field is 

H = Y, [j&i ■ Si+i - (-l)'D • Si x S i+1 

i 

- MS H ■ (g u + (-l)V) ■ Si] , (1) 

where the three terms are respectively the antiferromag- 
netic Heisenberg interaction, the DM interaction, and the 



Zeeman splitting energy. In addition to the superex- 
change coupling J > 0, D is the DM vector, H the 
external field, and g u,s are the uniform and staggered 
components of the gyromagnetic tensor. By making a lo- 
cal gauge transformation, which rotates the spins on two 
separate sublattices, and by neglecting all contributions 
at higher orders in D/J, Eq. ^ can be mapped to the 
simplified Hamiltonian 



H = 



JS, • Si+i - HSf - (-iyh s s? 



(2) 



where h s is an effective staggered field proportional to 
the product of H with a linear combination of D and g s \ 
novel features therefore arise only when such a material 
is subject to an external magnetic field. The Hamilto- 
nian of Eq. (2) can be mapped (by neglecting Zeeman 
splitting in the ^-direction) into a sine-Gordon model, 
a minimal framework whose bosonized version provides 
a good description of the opening of the spin gap^— 
Further unconventional features of this model include 
the specific heat£, magnetic susceptibility, 3,9 dynamical 
structure factor jSiLS line shape in electron spin resonance 
(ESR) measurements)* 4 magnetization^^ and the pres- 
ence of mid- gap states Air— 

While these numerous studies assumed that the sec- 
ond term in Eq. (2) has no significant role for weak mag- 
netic fields, in fact the uniform and staggered compo- 
nents may compete. A complete description of systems 
with DM interactions in arbitrary fields still requires 
the full exploration of the original Hamiltonian (1). A 
systematic investigation by Wang and coworkers 1 ^ using 
the density-matrix renormalization-group (DMRG) tech- 
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nique showed that the low-energy, high-field properties of 
copper benzoate are dominated by the uniform magnetic 
field, on which the spin gap depends linearly, while a 
crossover regime exists at intermediate fields. This pre- 
diction was later confirmed by ESR experiments in ap- 
plied fields up to 35 TM 

Most staggered-field studies have focused on ma- 
terials which are almost ideally ID in nature, such 
as copper benzoatefiiiSr— copper pyrimidine ^ 10 ' 19 " — 
Yb4As3j22r— and KCuGaFg. 25 However, real materials 
always have some interchain coupling, which in some 
cases may be comparable to the staggered magnetic field, 
if not also to the intrachain coupling. How the interchain 
interaction may change the essential physical properties 
remains an open question. Recent experiments on the 
weakly coupled chain system CuCi2-2(dimethylsulfoxide) 
(CDC) 26 ' 27 indicate that a gap opens at a finite value of 
h s , rather than at zero, when the uniform magnetic field 
is applied; the power-law dependence of the excitation 
gap seems to be different from that observed in quasi- 
1D materials. Early attempts to understand this behav- 
ior include calculations for the two-leg ladder; 28 ' 29 which 
show a QPT taking place as a consequence of the compe- 
tition between the staggered magnetic field and the inter- 
chain coupling. A chain mean-field theory developed^- to 
study the spin gap as function of the staggered field in 
2D and 3D found a spin gap opening immediately with 
the applied field (the "noncompetitive case" defined in 
Sec. II and discussed in Sec. IV), a conclusion confirmed 
by DMRG^i 

In a system where the staggered field competes with 
the ordering pattern favored by the magnetic interac- 
tions, a finite value of h s is required to induce a QPT. 
To date the only results available for this case are at the 
mean-field level, and may not deliver reliable conclusions 
for 2D systems. In particular, the linear dependence of 
the excitation gaps on the magnetic field appears to be in- 
consistent with experiment 2^ In this paper we contribute 
to the understanding of coupled Heisenberg chains in 
staggered magnetic fields by performing continuous-time 
Monte Carlo simulations using the worm algorithm^ - — 
and use our results to discuss the influence of interchain 
interactions and staggered fields on the ground state and 
the low-energy excitations. 

This paper is organized as follows. In Sec. II we pro- 
vide a formal introduction to the model and to the Monte 
Carlo method. In Sec. Ill we test our numerical tech- 
niques by computing the staggered magnetization in zero 
field, making contact with known results for coupled- 
chain systems and the square lattice. In Sec. IV we 
present our Monte Carlo results for the noncompetitive 
case, which demonstrate the 2D nature of the system. 
Section V contains our complete results for the compet- 
itive case, including the determination of the QPT, the 
staggered magnetization in the ordered phase, the gap 
in the disordered phase with corresponding fitting expo- 
nents, and a comparison of the phase diagram with differ- 
ent mean-field theories. We summarize our investigation 



in Sec. VI. 



II. MODEL AND METHOD 

In this study we investigate the S = 1/2 Heisenberg 
model on a spatially anisotropic square lattice in a stag- 
gered magnetic field. We stress that the purpose of our 
analysis is to determine, both qualitatively and quanti- 
tatively, the effects of a staggered field in different ge- 
ometries. We do not aim to make a direct comparison 
with experiment, but rather to demonstrate the phys- 
ical properties which may motivate the search for and 
characterization of an appropriate material in this class. 
Thus we focus here on effective models without the uni- 
form magnetic field. For completeness we consider both 
the "noncompetitive" case, where the geometry of the 
magnetic interactions and the staggered field are com- 
mensurate, and the "competitive" case, where they are 
not. We consider an antiferromagnetic intrachain cou- 
pling with both ferromagnetic and antiferromagnetic in- 
terchain couplings, which we will label respectively as 
{it, 0) and (ir, ir). The two different spatial arrangements 
of the staggered field, also (ir, 0) and (n, ir), then give one 
competitive and one noncompetitive situation for each 
case. 

The Hamiltonians can be expressed in the forms 

L L 
L 

+J2(-iy+ih s si P (3) 

i=l 

L L 

Hi 

i=l i=l 
L 

+ y £(- 1 )* h 'Stj, (4) 

i=l 

where L is the linear dimension of the lattice and J is 
the intrachain coupling, which we take as the unit of 
energy. The interchain coupling is Jj_, h s is the magni- 
tude of the effective staggered magnetic field, and Sij 
is the spin operator at lattice site with i the in- 

dex in the chain direction (x) and j the interchain in- 
dex (y-direction). The four situations are represented 
schematically in Fig. [TJ where the sign of Jj_ determines 
the nature of the interchain interaction and the dots or 
crosses correspond to up and down orientations of the 
applied staggered field. Figures [Ha) and (b) represent 
respectively the non-competitive and competitive cases 
of model ([3]), while Figs.QJc) and (d) represent the non- 
competitive and competitive cases of model dU). 

For our investigation of the ground states and low- 
est excitations of these models, we use the continuous- 
time worldline quantum Monte Carlo (QMC) method 
with "worm" updates. The concept of continuous-time 
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FIG. 1: Representation of couplings and staggered fields for 
the four possible cases considered here. Dots denote a field 
positive along z and crosses a negative field. The magnetic 
interactions have (n,ir) geometry in cases (a) and (d), and 
(it, 0) geometry in cases (b) and (c). The applied field ge- 
ometry is (n,n) in cases (a) and (b), which are described by 
Eq. $3$, and (tv, 0) in cases (c) and (d) [Eq. |@J]. Cases (a) 
and (c) are therefore non-competitive, while (b) and (d) are 
competitive. 



worldlines was first induced in QMC by N. V. Prokof'ev 
and coworkers in 1996/ 3 - Unlike the standard, discrete- 
time QMC algorithms, which are based on the Suzuki- 
Trotter decomposition, in the continuous-time method 
the Hamiltonian H = Hq + V is separated into a diago- 
nal term and a perturbation term. For a time-dependent 
perturbation, the partition function Z = Tr(e~ @ H ) can 

be expressed as Z = Tr[e~^°Z7(/3)], with U((3) the Mat- 
subara evolution operator, 



U(j3) = l - / dnV(n)U( n ) 



(5) 
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in which V(r) = e TH °Ve~ TH ° . In this series, each in- 
tegral corresponds to a worldline configuration with n 
"kinks," located at the points < r„ < r„_i < . . . < 
T\ < ft and varying continuously in the imaginary- 
time direction^ The integrals are evaluated by Monte 



Carlo sampling of the kink configurations. In contrast 
to discrete-time QMC algorithms, there is no systematic 
error caused by imaginary-time discretization. 

Based on this continuous-time worldline formulation, 
Prokof'ev and coauthors also developed an update algo- 
rithm based on breaking a worldline by inserting a pair 
of creation and annihilation operators, which then evolve 
by local moves^ (This became known as the "worm" up- 
date from the motion of pairs of worldline kinks.) The 
algorithm considers two configuration spaces, denoted Z 
and G: while Z contains only closed worldlines, and is the 
space of the partition function, G contains a worldline 
broken by a physical process connecting points (fi,Ti) 
and (fj, Tj). The Z and G configuration spaces can be ex- 
changed by creation or annihilation of an (i,j) "worm" 
pair on the same fiat worldline. This is equivalent to 
processes in the grand canonical ensemble. In the G con- 
figuration space, processes which move the worms in both 
time (changing worm position in the imaginary time di- 
rection on the same site) and space (the worm is moved 
to a neighboring site, creating or annihilating a kink) 
update the configuration and change its structure. 

In the Z configuration space it is easy to compute 
thermodynamic quantities such as energy, magnetization, 
susceptibility, and specific heat. In the G configuration 
space, each accepted worm move results in a contribution 



to the histogram of the Green function, G(i 



■Tj). 



With sufficiently many Monte Carlo steps, one obtains a 
convergent Matsubara Green function, defined as 

G(r, r) = G(fi - r 3 ,n - r 3 ) = (T T S+(n)SZ ( Tj )), (6) 

where T T is is the time-ordering operator. The average 
(. . .) is obtained from the histogram, and is normalized 
by G(0,0). From the statistical average of the Green 
function, one may obtain the critical exponents of G, the 
single-particle excitation spectrum, and certain correla- 
tion functions. 

In comparison with previous results by the DMRG 
technique, 31 QMC algorithms have the major advantage 
of treating 2D models. In the continuous-time worldline 
approach, lattice sizes of L up to 100 and inverse tem- 
peratures (3 = 1/T up to 100 can be accommodated in 
approximately 24 hours of calculation time when working 
on a server with an Inter XEON E5460 CPU. A further 
advantage of QMC is that it is easy to work with peri- 
odic boundary conditions, so that there are no difficulties 
caused by edge states^ We therefore expect that reliable 
numerical results can be obtained for the 2D case by this 
technique. In this paper, we focus on the energy gap A, 
which characterizes gapped phases, and the transverse 
staggered magnetization M£, which is characteristic of 
ordered phases. The gap is obtained 3 ^ from the Green 
function through the expression 

A _ ln[G(p,T)/G(p,r )} 
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while the staggered correlation function 
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FIG. 2: (Color online) Staggered magnetization Mj_ as a 
function of interchain coupling for both signs of Jx, at zero 
staggered field. The isotropic square lattice is re prese nted 
by Jx = 1. Inset: Mj_ as a function of — yj Jx|(l + 
0.095|J ± |)ln 1/3 |J ± /1.3|,— showing the strong logarithmic 
corrections to the square-root form of a chain mean-field the- 
ory. Up- and down-pointing triangles denote data obtained 
for LxL samples with L = 100, while squares and diamonds 
denote data obtained on rectangular samples with aspect ra- 
tio Lx/L v = 8. The solid line is the fitting function deduced 
in Ref. M. 



is obtained from the Green-function histogram. 
The staggered magnetization is defined as 

(tt, 7r) for model ((3|), 



and Q 



ErQjM/iV, where Q 
(tt, 0) for model (JJ). 



III. MAGNETIZATION AT ZERO FIELD 

We begin by investigating the case h s — 0, which 
also serves to benchmark our method for accuracy 
and systematic errors. We consider both the spatially 
anisotropic antiferromagnetic Heisenberg square lattice 
and the lattice of Heisenberg chains with ferromagnetic 
interchain coupling. For \ J±\/J > 0.15, the lattice size L 
in our calculation is 100, with j3 set equal to L. For lower 
values of the interchain coupling, we find that rectangu- 
lar lattices are required to achieve reliable results, and we 
adjust the cluster aspect ratio accordingly. The staggered 
magnetization we compute for the 2D antiferromagnetic 
Heisenberg model is shown in Fig. [2] For benchmarking 
purposes, in the isotropic square lattice (J = J± = 1) 
we find Mj_ = 0.313 ± 0.006, which agrees very well with 
the known value of 0.306.i°- The case of ferromagnetic 
interchain coupling is not entirely symmetrical, the non- 
universal behavior at larger values of | J± | being manifest 
as a more rapid saturation in the antiferromagnetic case. 

The problem of the weakly coupled Heisenberg spin 
chains is a fundamental one in quantum magnetism, and 
has received a great deal of attention over the last five 



decades. It encapsulates the physics of the crossover 
from truly ID systems, dominated by quantum fluctu- 
ation effects, to high-dimensional, renormalized classical 
behavior. Because the S = 1/2 Heisenberg chain is al- 
ready critical, with a gapless ground state, any trans- 
verse coupling in an unfrustrated geometry is thought 
to give rise to magnetic order, and extensive investiga- 
tion has reinforced the general agreement that the critical 
J± c = 0. Early discussion of the critical behavior of the 
ordered moment and Neel temperature used the high- 
dimensional framework of renormalized spin waves, and 
suggested a purely logarithmic rise of M± with |Jj_|^ 
By contrast, in a chain-based weak-coupling approach, 36 
where the transverse interactions are modeled as an ef- 
fective staggered magnetic field, the ordered moment has 
a power-law dependence, M± oc \J\ J±\, albeit with sus- 
pected logarithmic corrections. The definitive study of 
the problem was performed by Sandvikj 37 using QMC 
within a multichain mean-field theory, and reveals strong 
logarithmic corrections to the square-root dependence, 
M± oc — V^[Jl| In 1 / 3 | Jx (with a weak additive linear 
term) . The power of 1 /3 in the logarithm was also found 
in the problem of the single chain in a staggered fields 
The importance of the logarithmic terms lies in the pres- 
ence of marginally irrelevant operators, which are ne- 
glected in the transformation of the high-dimensional 
Heisenberg model to a chain in an effective staggered 
field. 

Our data for both signs of the transverse coupling, 
shown in Fig. [21 have the same form at small |Jj_|. Our 
results are fully consistent with those of Ref. [33]. First, 
we confirm our sensitivity to the expected logarithmic 
corrections, which will be essential in the sections to fol- 
low. Second, our results for square lattices deviate from 
the expected form at values of the transverse coupling 
below | Jj_| « 0.15. We simulate instead rectangular sys- 
tems of different aspect ratios up to L x /L y = 8,^ 7 which 
allows us to obtain accurate values of Mj_ at least to 
| J±\/J = 0.05, below which the calculations become very 
time-consuming. This study of the staggered magnetiza- 
tion at h s = therefore allows us to benchmark the ac- 
curacy of our results for highly anisotropic systems. The 
effect of the finite temperature in our calculations is not 
thought to be significant. 37 Quantifying the logarithmic 
corrections to the chain mean-field picture in this way 
will be important in Sec. V, where we will use it as a 
guide to understanding our numerical results in the pres- 
ence of a staggered field. We will show that the field acts 
to alter significantly the effective dimensionality of the 
system and the relevance of these logarithmic terms. 



IV. NONCOMPETITIVE CASE 

In the two noncompetitive cases, the magnetic inter- 
actions and the applied staggered field have the same 
spatial pattern [Figs. 1(a) and (c)]. For any finite value 
of Jj_, the system is magnetically ordered at h s = 0, by 
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FIG. 3: (Color online) Excitation gap as a function of stag- 
gered magnetic field for the 2D spin-i Heisenberg square lat- 
tice (J± = J) in the noncompetitive case. Monte Carlo results 
are presented as filled circles for the (7T, 7t) geometry [model 
@] and as squares for (71", 0) [model @]. The solid lines are 
fitting curves. Shown also (triangles) are Monte Carlo results 
for the Heisenberg chain. 



FIG. 4: (Color online) Excitation gap as a function of stag- 
gered field for the anisotropic 2D spin-i antiferromagnetic 
Heisenberg model in the noncompetitive case for J±/ J — 0.1. 
Triangles and squares correspond respectively to h 3 steps of 
0.005 and 0.2. The fitting curve is obtained from the data 
with Sh s — 0.005 (inset). The blue, dashed line marks the 2D 
regime for J±/J = 0.1. 



a spontaneous breaking of the SU(2) spin symmetry, and 
has massless spin- wave excitations. When a commensu- 
rate staggered field is applied, the symmetry is broken 
explicitly as the field direction selects the spin orienta- 
tion. As for a ferromagnet in a uniform field, the stag- 
gered field serves as an anisotropy term, which opens a 
spin- wave gap for any finite value of h s . In the conven- 
tional mean-field approach^ a standard linear spin- wave 
approximation gives a gap 

A « y/4SJh s [l + (h s /ASJ)]. (9) 

By extrapolation from smooth and rapidly convergent 
DMRG results for antiferromagnetic Heisenberg ladders 
in staggered fields, the authors of Ref. [3l| also obtained 
a 2D field-dependence of A = (2.27 ± 0.01) ft0.50io.oi_ 
However, the noncompetitive geometry may also be dis- 
cussed within a chain mean-field theory, where the inter- 
chain interactions are treated as an effective staggered 
field, h = —2J±M±(h)r^ which is reinforced by the ap- 
plied staggered field. In this treatment, the dominant 
physics is the breaking of the continuous symmetry of 
the spin direction, and the exact dependence on J± is 
not clear; that the gap A is independent of J± in Eq. flS]) 
underlines the fact that it is in essence purely an effect 
of the field on the chain, which would suggest some in- 
fluence of ID physics. 

Our analysis is the first direct numerical calculation of 
the gap evolution in this case. We begin by considering 
the "isotropic" square lattice (| J±\ — J) with systems of 
L = 100 and temperatures f3 = 100, to maximize the reli- 
ability of our data by systematic extrapolation. In Fig. [3] 
we show the dependence of the gap on the staggered field 
for both noncompetitive cases. A fit to the formula 

A = a h* (10) 



yields excellent results (Fig.©, with A = 2.19ft°- 50±a01 
for the (n,ir) case [Fig. 1(a)] and A = 1.65ft°- 50±0 Q1 for 
the (7r,0) case [Fig. 1(c)]. The exponents in both cases 
agree perfectly with the result predicted both by mean- 
field theory and by extraolation from Heisenberg lad- 
ders. The isotropic square lattice shows unambiguously 
2D behavior, with an immediate opening and square-root 
growth of the gap. 

By contrast, in the purely ID case is it knowr*^^ that 

2/3 . 

the gap opens according to A cx h s with logarithmic 
corrections. This behavior is shown as the triangles in 
Fig. [3J The more rapid growth of the gap in 2D may 
be considered heuristically as the consequence of a mu- 
tual reinforcement of the applied and effective staggered 
fields, the latter arising from the magnetic moment en- 
hanced by the former. The key question to address 
is whether the system displays any kind of continuous 
crossover, as a function of |J_l|, from a regime charac- 
terized by ID exponents to a 2D regime. By taking the 
non-competitive (tt, it) case at J± = 0.1 and fitting the 



TABLE I: Fitting coefficients and exponents for five data sets 
with different point spacings Sh s . 
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FIG. 5: Schematic representation of the spin state in the com- 
petitive case on varying the staggered field h 3 . The system 
has FM interchain interactions and a (it, it) staggered field 
[model ©]. Black arrows represent the spin directions and 
the staggered field is represented by dots (up) and crossed 
circles (down). 



gap obtained from ten sets of data points with five dif- 
ferent h s intervals, 0.005, 0.01, 0.05, 0.1, and 0.2, we find 
the fitting parameters listed in Table [J] It is clear that 
the exponent changes from a 2D form to a ID form: for 
large h s , the system is effectively no longer aware of the 
coupling, and has ID behavior, while at sufficiently small 
h s the 2D behavior always emerges. The same result is il- 
lustrated in Fig. fusing the data points with 5h s = 0.005 
and 5h s = 0.2. It is clear that the latter data set does 
not fall on the square-root curve obtained from the for- 
mer. We may conclude that the 2D regime in this model 
is given approximately by h s < Jj_/2. 



V. COMPETITIVE CASE 

We turn now to the competitive case, where the in- 
terchain coupling and the staggered field compete to es- 
tablish the pattern of magnetic order. From a mean- 
field analysis of this system^ there exist two different 
phases at finite h s , an ordered phase with spontaneous 
symmetry-breaking (SSB, of the continuous U(l) sym- 
metry) in the plane normal to the staggered field and a 
gapped, "symmetric" phase in which the spins are ori- 
ented in the field direction. 

A phenomenological description of the situation is il- 
lustrated in Fig. [5] for the case of a (it, n) staggered field 
applied to a system with J± < 0. At h s = 0, the system 
is ordered with a spontaneous breaking of SU(2) sym- 
metry and one-site translational invariance in the trans- 
verse direction. For h s — + , a spin-flop transition occurs 
into the plane perpendicular to the field, represented as 
the (xy) plane, but there remains a SSB in this plane 
[Fig. Ufa)]. For finite h s , the spins are forced to rotate 
into the direction of the staggered field, or into the yz 
plane in Fig. [5^b) , adopting a canted structure with two- 
site translational invariance in the transverse direction. 
The finite magnetic order parameter in this phase is re- 
duced by the staggered field, and the excitations of the 
system remain gapless. Finally, when the staggered field 
is increased beyond a critical value h c , the order param- 
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FIG. 6: (Color online) Staggered correlation function C s (r) = 
(S^S~e lQ - T ) (see text) with L = 64 and J ± /J = -1.0, for 
two staggered fields, h s = 1.05 and 1.2, chosen to represent 
respectively the ordered and the gapped phase. 



eter is driven to zero, long-range order collapses, and the 
spins are fully oriented in the field direction [Fig.[5][c)]. 

We begin by our numerical analysis by illustrating the 
staggered correlation function within the system. Figure 
[6] shows C s (r) as a function of distance for model (0]) , 
with J± = — 1 and for two values of h s . For h s = 1.05, 
staggered correlations are finite and there is long-range 
order in the plane transverse to the applied field. By con- 
trast, for h s — 1.2, C s (r) falls to zero abruptly away from 
the edges of the system; these two values of h s therefore 
fall on opposite sides of the anticipated critical staggered 
field, h c . We use the staggered magnetization to analyze 
both the SSB phase and the field-induced quantum phase 
transition between the SSB and symmetric phases. We 
characterize the symmetric phase by its gap, defined in 
Eq. ([7]), and by the scaling of this gap. 



A. Critical point 

From the discussion above, the critical point h c is the 
single most important quantity in the descriptiontion of 
the competitive case. Once h c is determined, a process 
achieved most accurately using the behavior of the order 
parameter in the SSB phase, the gap and the 2D excita- 
tions of the symmetric phase can be calculated with high 
accuracy. To determine the critical point, or the phase 
boundary in the space of J± and h s , we apply the finitc- 
size-scaling hypothesis to the behavior of the staggered 
correlation functions C^o and C T .„ [respectively for mod- 
els and (H])] as functions of L. These are expected^ 
to obey the scaling form 



C(8) = L 2 - d - z ](bL x l v \P/L Z ), 



(11) 



where 5 = h s — h Cl f3 is the inverse temperature, and z is 
the dynamical exponent. 

The critical point h c can be measured accurately by 
computing the staggered correlation function near the 
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FIG. 7: (Color online) Determination of h c from the finite- 
size-scaling hypothesis for (a) J± = 0.5 and (b) J± — —0.5. 
Error bars are set by the inexact crossing of the lines, shown 
in the insets. In panel (a), 0.6471 < h s < 0.6472 and 0.462 < 
LCr^ < 0.471. In panel (b), 0.49300 < h s < 0.49315 and 
0.502 < LGx , < 0.510. 



critical point for different lattice sizes L x L. On fixing 
f3/L z = 1, the scaling function / depends on a single 
parameter, SL 1 ^. Precisely at the critical point, <5 = 0, 
CL d+z ~ 2 is independent of the lattice size, which implies 
that calculated curves for CL d+z ~ 2 as functions of h s for 
different lattice sizes should cross at the critical point. It 
is clear from this form that z is a measure of the system 
coherence: as z becomes larger, the coherence of the sys- 
tem vanishes more rapidly in space at the critical point. 
Here we test that the dynamical exponent takes the value 
z = 1, and we find this to hold at all points around the 
critical region. 

We have simulated 11 different values of | J±] for both 
ferro- and antiferromagnetic interchain coupling, with 
fixed L/f3 = 1 and at least four lattice sizes, namely 
L = 40, 48, 60, and 64. We have verified the accuracy 
of our results by performing calculations with L = 100 
for the case J± = 0.5. For a determination of the critical 
point, we work only with square samples at all values of 
JjJ; the simultaneous vanishing of the superfluid den- 
sity (spin stffness) in both directions at h c negates the 
advantanges of rectangular samples for low | J±\ (Sec. Ill) 
in this case. Examples of the sets of crossing curves are 
illustrated in Fig. [7] for J± = 0.5 and —0.5. The error in 
h c can be estimated from the separation of the different 
intersection points within a single manifold of curves; all 
errors Sh c are of order 10 -4 . The full numerical details 
determined from this procedure are presented in Table 
ITIl The phase diagram deduced from these values of h s 
is discussed in Sec. VD. 



B. Staggered Magnetization 

For applied fields h s < h c , the SSB phase possesses a 
transverse staggered moment (or off-diagonal long-range 
order) Mj\ We have calculated this quantity for the 
case h s = in Sec. III. In general one expects the mag- 
netic order to be suppressed by the competing staggered 
field, because of the different ordering patterns they fa- 
vor, until Mj_ vanishes at h s = h c . In Fig. [5] we show 
Mj_ as a function of h s for different values of the in- 
terchain coupling J±. For larger values of J±, of both 
signs, the order parameter shows a conventional, mono- 
tonic decrease [Fig.[5Ja)]. However, when |«7j_ | < 0.2, for 
both signs of J± [Figs. E^b) and (c)], we observe that 
first increases despite the increase in the compet- 
ing staggered field. We stress that our results for small 
\J±\ were obtained on rectangular samples with aspect 
ratio L x / L y = 8, following the procedure established in 
Sec. Ill as providing the highest available accuracy. 

To our knowledge, this novel and purely quantum 
mechanical effect has not been remarked upon previ- 
ously in coupled Heisenberg chains. An analogous ef- 
fect can be found in dimerized quantum magnets 4 ^ such 
as NH4CUCI3, where an applied uniform magnetic field, 
which in principle competes with the transverse staggered 
order, nonetheless causes the staggered moment to rise at 
the same time as the uniform polarization is increased. 
A heuristic understanding of this phenomenon is that 
all forms of magnetic order are in fact competing with 
quantum fluctuation effects favoring disordered states. 
When the disordered state is suppressed by an applied 
field, more ordered spin "weight" is available both for 
the magnetization component favored by the field and 
for the component favored by the magnetic interactions. 

As h s is increased further, Mf is suppressed for all 



TABLE II: Critical staggered fields h c for systems with trans- 
verse couplings I J± \ varying from 0.05 to f. Shown are results 
for both ferromagnetic interchain couling in a (n, ir) staggered 
field ((3| and antiferromagnetic interchain coupling in a (n, 0) 
field. 



(tt,tv) 


(tt,0) 


J± 


h c 


5h c 


J± 


h c 


Sh c 


-0.05 


0.02459 


0.00007 


0.05 


0.02606 


0.00008 


-0.1 


0.06336 


0.00011 


0.1 


0.06964 


0.00027 


-0.2 


0.15805 


0.00015 


0.2 


0.18495 


0.00013 


-0.3 


0.26495 


0.00016 


0.3 


0.32470 


0.00022 


-0.4 


0.37775 


0.00014 


0.4 


0.48034 


0.00008 


-0.5 


0.49305 


0.00012 


0.5 


0.64715 


0.00015 


-0.6 


0.60990 


0.00006 


0.6 


0.82166 


0.00018 


-0.7 


0.72630 


0.00012 


0.7 


1.00225 


0.00020 


-0.8 


0.84300 


0.00011 


0.8 


1.18617 


0.00014 


-0.9 


0.95850 


0.00017 


0.9 


1.37339 


0.00015 


-1.0 


1.07300 


0.00021 


1.0 


1.56290 


0.00019 
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FIG. 8: (Color online) Transverse staggered order parameter 
M± as a function of h s for different values of the interchain 
coupling Jj_. (a) Regime of large |Jj_|, showing conventional 
suppression of Mj_ by the competing field h s . (b) and (c) 
Regime of small \J±\, showing unconventional increase of M ± 
with h s at low applied fields. 



values of \J±\, and it falls continuously to zero at the 
second-order quantum phase transition. In this regime 
the finite-size effects in our simulations are large, even 
though the statistical errors are very small. Unlike the 
determination of h c itself, there is neither a systematic 
approach by extrapolation with different system sizes, 
which may allow a sufficiently accurate determination of 
M x close to h c , nor a method for the accurate estima- 
tion of such errors. Thus we are not able to deduce the 
critical exponents of the transverse staggered moment in 
order to characterize this side of the quantum phase tran- 
sition. The results of Sec. VC below suggest that some 
anomalous behavior may be expected. 



C. Energy gap 

In the gapped, symmetric phase we wish to consider 
the properties of the excited states. We compute the 
energy gap A and extract the scaling behavior of the 
gap in the region h s ~h c < \ J±\/2 deduced in Sec. IV. In 
Fig. Owe present selected examples of the function A(h s ) 
for a range of positive and negative values of J± . Unlike 
in the SSB phase, there are no strong finite-size effects 
on the gap data close to h c for the symmetric phase, a 
result we have confirmed by studies on lattices of various 
sizes up to L = 100. For a given /i c , determined from 
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FIG. 9: (Color online) Gaps as a function of staggered mag- 
netic field h s for four different values of | Jx \ ■ Monte Carlo re- 
sults for the gap are presented by black circles and red squares 
respectively for ferro- and antiferromagnetic interchain cou- 
plings. The critical points h c are determined from finite-size 
scaling of the transverse order parameter in the SSB phase. 
The gap is fitted to the form A = ao(h a —h c ) a . Fitting curves 
are shown by solid lines and complete fitting parameters pre- 
sented in Table mil 



finite-size scaling of the transverse order parameter in 
the SSB phase, we collect eight data points within the 
scaling regime h s — h c < \J±\/2. We fit our numerical 
data to the formula A = ao(h s — h c ) a . 

Our gap calculations are of course consistent with the 
values of h c computed in the SSB phase. The gaps for Jj_ 
values of opposite sign show only small, quantitative dif- 
ferences. The two gaps converge to that of the ID chain 
when J± — > 0. As in Sec. VA, h c increases with the mag- 



TABLE III: Fitting coefficients and exponents for the excita- 
tion gap in systems with traverse couplings | J± | varying from 
0.05 to 1. Shown are results for both ferromagnetic interchain 
coupling in a (ir, 7r) staggered field © and antiferromagnetic 
interchain coupling in a (tt, 0) field Q. 



(7I",7r) 


O,o) 


J± 


ao 


a 


5a 


Jx 


do 


a 


5a 


-0.05 


0.929 


0.693 


0.003 


0.05 


0.943 


0.698 


0.004 


-0.1 


0.926 


0.718 


0.005 


0.1 


0.880 


0.708 


0.010 


-0.2 


0.862 


0.724 


0.006 


0.2 


0.822 


0.715 


0.005 


-0.3 


0.815 


0.718 


0.006 


0.3 


0.759 


0.706 


0.008 


-0.4 


0.761 


0.704 


0.004 


0.4 


0.744 


0.703 


0.003 


-0.5 


0.735 


0.700 


0.004 


0.5 


0.732 


0.702 


0.006 


-0.6 


0.716 


0.693 


0.003 


0.6 


0.686 


0.692 


0.008 


-0.7 


0.709 


0.694 


0.005 


0.7 


0.676 


0.689 


0.009 


-0.8 


0.705 


0.692 


0.005 


0.8 


0.673 


0.688 


0.005 


-0.9 


0.698 


0.692 


0.007 


0.9 


0.667 


0.686 


0.006 


-1.0 


0.693 


0.692 


0.007 


1.0 


0.661 


0.686 


0.008 
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FIG. 10: (Color online) Critical exponents a of the gap A 
(Table UTT)) as a function of \Jj_\ for (a) ferromagnetic inter- 
chain coupling and (b) antiferromagnetic interchain coupling. 
The solid lines are a guide to the eye. 



nitudc of J± , reflecting the fact that it is the competition 
between the staggered field and the interchain coupling 
that drives the quantum phase transition. Full details of 
the fitting parameters for the gap, including the errors in 
the exponents we extract, are listed in Table ITTT1 

Figure [TU] shows the gap exponent extracted for differ- 
ent values of J±. Because all of the data is fitted to a 
simple power law, the logarithmic corrections in the ID 
chain result in an effective exponent of 0.61 for J± = 0. 
The exponent shows an initial increase as | Jj_| increases, 
but then falls weakly and remains close to a — 0.7 over 
the remainder of the range up to \J±\ = 1. The values 
for the two cases of ferro- and antiferromagnetic J± dif- 
fer only very slightly, and no universal value is indicated. 
It is clear, however, that the value of a does not fall to 
0.5, even arbitraily close to the transition, as would be 
expected in a 2D mean-field theory or from the results 
of Sec. IV. Our results therefore provide strong evidence 
for nontrivial quantum physics. 

For a heuristic explanation of these data, we appeal to 
the chain mean- field theory of Ref. (3(| , where the inter- 
chain interactions are treated as an effective staggered 
field ho = — 2J±M ( [. The presence of a real staggered 
field either reinforces this effective one, as in the non- 
competitive geometry of Sec. IV, or suppresses it in the 
competitive geometry Unlike the noncompetitive case, 
where the applied field breaks all continuous symmetries, 
in the competitive case a U(l) symmetry is maintained 
and the system may still be treated as a chain in a single 
effective field. 

At lowest order, one may write h e e = —2J±M±(h s ) = 
-2J±M^ + h s , from which M±(h s ) = M^(l - 
h s /2J±M c [). This simple linear approximation contains 
directly the competition between the applied and effec- 
tive staggered fields, and suggests a quantum phase tran- 
sition at h c = 2J±Mj_, where the applied field cancels the 
interactions. The quasi-linear relation between h c and 
J± is clear from our calculations, while the relevance of 
next-order, and possibly also logarithmic, corrections is 
clear from our results with ferro- and antiferromagnetic 
J± . The spin chain in a positive effective staggered field 



is once again a problem with no continuous symmetry, 
in which a spin gap opens directly. The key point about 
this heuristic picture is that the cancellation of the ap- 
plied and effective staggered fields results in a quasi-lD 
problem close to h c , and hence the anomalous exponents 
we observe for all values of J± are to be expected: h c 
is always the staggered field required to cancel the 2D 
coupling and reduce the model to a ID system. 

The inexact values of the exponent can be ascribed 
to departures from universality in the ID nature of this 
system, and are are not indicative of a phase genuinely 
intermediate between ID and 2D. Because of the effec- 
tively ID nature of the system close to the critical applied 
staggered field, logarithmic corrections to the gap scaling 
may remain significant for all values of | J± \ . We conclude 
these considerations by commenting on the coupled spin- 
chain system CDC,— where anomalous gap exponents 
have been measured despite the fact that the interchain 
coupling is thought to be significant. We suggest on the 
basis of our results that the origin of this behavior may 
lie in the effects of a competitive staggered field. 



D. Phase diagram 

Finally, we present in Fig.lTTlthe phase diagrams of the 
two competitive-case models for ferro- and antiferromag- 
netic interchain coupling. The data is taken from Table 
II. For comparison we show also the chain mean-field re- 
sults and 2D mean-field results of Ref. [3(| ■ The numeri- 
cally exact boundary determined by our quantum Monte 
Carlo simulations lies between those obtained from the 
two mean-field theories [Fig. fTTT b)]. showing that nei- 
ther is particularly accurate for the problem of the 2D 
staggered magnetic field, and confirming the general de- 
parture from universality of this type of system. 

The phase boundary, or critical point as a function of 
the interchain coupling constant, can be fitted well over 
the full range of J± by a simple power law, 



h c = b \J±\ 



(12) 



with parameters h c = 1.08|Ji| 1,15 for the case of ferro- 
magnetic interchain coupling and h c = 1.57|Jj_| for 
antiferromagnetic coupling. 



VI. SUMMARY 

In summary, we have studied the zero-temperature 
phase diagram and low-energy spin excitations of 
anisotropic, two-dimensional spin-^ Heisenberg models 
on the square lattice under a staggered magnetic field. 
We have used a continuous-time Monte Carlo method 
to calculate ordered moments and spin gaps, and hence 
phase boundaries and scaling relations. At zero field, 
we compute the properties of the anisotropic Heisenberg 
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FIG. 11: (Color online) Phase diagrams for competitive-case 
models in (a) (n, n) and (b) (n, 0) staggered fields, i.e. with (a) 
ferromagnetic and (b) antiferromagnetic interchain coupling. 
Monte Carlo results, given in Table (TT] are shown by solid 
circles and power-law fits (see text) by solid lines. In (b), the 
blue, dotted line marks the phase boundary predicted by a 
2D mean-field theory and the red, dashed line that from a 
chain-based mean-field theory.— 

model and benchmark our calculations for the complex- 
ities of high anisotropy and strong logarithmic correc- 
tions. When the applied field and magnetic interactions 
cooperate, a gap opens immediately and scales exactly as 



the square root of the staggered magnetic field. When the 
two compete, we find a field-driven quantum phase tran- 
sition from a gapless phase of staggered magnetic order, 
and characterized by spontaneous symmetry-breaking in 
the transverse direction, to a gapped, disordered, field- 
dominated phase. In the ordered phase, we find a novel 
enhancement of the staggered moment by a competing 
field, a purely quantum effect. We determine scaling 
regimes, discuss the scaling properties of the gap, and 
show that the physics of the system in the competitive 
case is essence always one-dimensional at the quantum 
critical point due to the cancellation of applied and ef- 
fective staggered fields. 
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